Vortex nucleation and quadrupole deformation of a rotating Bose-Einstein condensate 
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Using a semi-analytic model based on the Thomas-Fermi approximation, we investigate the 
relevance of the quadrupole deformation of a trapped Bose-Einstein condensate for the nucleation of 
quantized vortices. For sufficiently high angular velocities SI of the trap, the tendency of the system 
to exhibit spontaneous deformation is shown to lower the barrier which inhibits the nucleation of 
vortices at smaller fl. The corresponding value of the critical angular velocity fl c is calculated as a 

' function of the deformation of the trap and of the chemical potential. The theoretical predictions 

for Q c refer to the case of a sudden switch-on of the deformed rotating trap and they are compared 

. with recent experimental data. 
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^ ■ I. INTRODUCTION 
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The problem of the nucleation of quantized vortices in Bose superfluids has been the object of recent experimental 
work with dilute gases in rotating traps It now emerges clearly that the mechanism of nucleation depends 

crucially on the actual shape of the trap as well as on how the rotation is switched on. A first approach consists 
of a sudden switch-on of the deformation and rotation of the external potential which generates a non-equilibrium 
t^J- • configuration Jl|-Q]. For sufficiently large values of the angular velocity of the rotating potential one observes the 
£SJ ' nucleation of one or more vortices. In a second approach either the rotation or the deformation of the trap are 
switched on slowly ||,|5| . In this way, the system is in conditions of equilibrium up to some critical velocity above 
which it exhibits a dynamical instability giving rise to the nucleation of vortices. In the case of a sudden switch-on 
of the deformation and rotation of the trap, the observed critical angular velocity turns out to be close to the value 
associated with the instability of the surface oscillations of the condensate. This occurs at the angular velocity |^,|7j 
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. where ui(£)/2ir is the frequency of the surface oscillation with angular momentum £ and we have assumed that the 
deformation of the trap is negligibly small. The actual value of I fixing the critical angular velocity depends explicitly 
on the shape of the rotating potential and, to some extent, can be selected in the experiment Eq. (Q) corresponds 
O • to the famous Landau criterion for superfluidity applied to the problem of rotations. In the Thomas- Fermi limit the 
critical angular velocity takes the simple form Q u> CI = loj_ / where uj± is the transverse frequency of the harmonic 
confinement. The experimental evidence that Eq. (|l|) provides a good estimate for the critical angular velocity for 
the nucleation of quantized vortices points out the crucial role played by the shape deformation of the condensate, 
in agreement with the theoretical considerations developed in fi|il3|-|l3|] . This is further confirmed by the direct 
observation Hl^] of strong deformations occurring temporarily during the process of nucleation. 

The purpose of this paper is to gain further insight into the mechanism underlying vortex nucleation when the 
deformation and rotation of the trap are switched on suddenly. So far, the explicit relation between the deformation 
instability and the vortex creation was not enough clear. In the absence of shape deformation the nucleation process 
is inhibited by the occurrence of a barrier located near the surface of the condensate. Using the Thomas-Fermi (TF) 
approximation to the Gross-Pitaevskii (GP) theory at zero temperature |l4|] we will show that this barrier is lowered 
by the explicit inclusion of shape deformations and eventually disappears at sufficiently high angular velocities, making 
it possible for the vortex to nucleate. We will consider the simple but relevant case of quadrupole deformations which 
has already been the object of systematic experimental studies [jl],||-f|. 

II. VORTICES IN AXISYMMETRIC CONFIGURATIONS 

A quantized vortex is characterized by the appearance of a velocity field associated with a non- vanishing, quantized 
circulation. If we assume that the vortex can be described by a straight line the quantization of circulation takes the 
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simple form 

Vxv vortcx = — <5^(r-d)i (2) 
m 

for a vortex located at distance d = |d| from the z-axis. The general solution of ^ can be written in the form 

Vyortex = V (<p d + S) (3) 

where is the azimuthal angle around the vortex line at position d and S is a single-valued function which gives 
rise to an irrotational component of the velocity field. The irrotational component may be important in the case of 
vortices displaced from the symmetry axis and its inclusion permits to optimize the energy cost associated with the 
vortex line |I5| ] . Considering a straight vortex line is a first important assumption that we introduce in our description 

B 

The inclusion of vorticity is accompanied by the appearance of angular momentum and by an energy cost. Under 
suitable conditions the system may nevertheless like to acquire the vortical configuration. This happens if there is a 
total energy gain in the rotating frame where the system is described by the Hamiltonian 

ff(fi) = H - flL z . (4) 

Here H is the Hamiltonian in the laboratory frame, L z is the angular momentum, and Q is the angular velocity of 
the trap around the z-axis. 

The simplest way to calculate the angular momentum associated with a displaced vortex is to assume axi-symmetric 
trapping and to work with the TF-approximation. In this limit the size of the vortex core is small compared to the 
radius of the condensate so that one can use the vortex-free expression n(r) = fi[l — (r±/R±) 2 — (z/Z) 2 ]/g for the 
density profile of the condensate, where r\ = x 2 + y 2 and g = 4irah 2 /m is the coupling constant, fixed by the positive 
scattering length a. Then, one can write the angular momentum in the form 

L z = m Jdz Jdr±r±_n(rx,z) ®v vortex • dl , (5) 

where the line integral is taken along a circle of radius r±. Use of Stokes' theorem gives the result M 



L z (d/R±) = Nh 
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where d is the distance of the vortex line from the symmetry axis. Eq. (^|) shows that the angular momentum per 
particle is reduced from the value h as soon as the vortex is displaced from the center. 

To calculate the energy cost associated with the vortex line one should in principle start out with Eq. (^) and 
optimize the choice for the phase S. In order to obtain a simplified description we will make use of the expression 



E v {d/R ± ,fi) = E v (d = 0,^) 
which generalizes the TF-result mM 
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E v (d = )f i) = — Z\og[ — ^ =7Vf^ ^ fog 1.342-^ . (8) 

for the energy of an axi-symmetric vortex. Here, Z is the TF-radius in z-direction and £o is the healing length 
calculated with the central density uq. The d-dependence contained in (Q) is simply understood by noting that the 

factor 4no^ [l — (d/R±) 2 ~\ 3 ^ 2 /3 corresponds to the column density Jdz n(d, z) evaluated with the TF-approximation. 
Expression (0) is expected to be correct within logarithmic accuracy (see also jl9| and references therein). It could 
be improved by including an explicit c?-dependence in the healing length inside the logarithm, in order to account for 
the density dependence of the size of the vortex core. 

Eqs. (Q) and (Q) show that when d = R± neither angular momentum nor excitation energy is carried by the system. 
Of course these estimates, being derived with the TF-approximation, are not accurate if we go too close to the border. 
If the surface of the condensate is described by a more realistic density profile both the angular momentum and the 
energy are nevertheless expected to vanish when the vortex line is sufficiently far outside the bulk region. Within the 
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simplifying assumptions made above we can conclude that the configuration with d = R± corresponds to the absence 
of vortices, while the transition from d = R± to d = describes the path nucleating the vortex. The energy associated 
with each point of this path in the rotating frame is simply calculated using the Hamiltonian (^). It is given by Jl7 19 



E v (d/R±,n,(x) = E v (d = 0,n) 
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Some interesting features emerge from Eq. . First one finds that the occurrence of a vortex at d = is energetically 
favorable for angular velocities satisfying the condition fl > £l v ([i) = E v (d = 0, fi)/NH. This is the well known criterion 
for the so called thermodynamic stability of the vortex. It is worth noticing that in the TF-limit one should have 
tt v (/i) /ujj_ << 1. In fact, using the result /i = griQ for the chemical potential and the expression /i = mu>'j_R 2 [_/2, one 
can write 



w± 2 \R_lJ V Co 



(10) 



which tends to zero when R± ^> a±. Here, a± — ^/h/miv± is the radial oscillator length. In the actual experiments 
the ratio ( |l0| ) is not very small. For example, in the case of Ref. M Q v ~ 0.35 u>±. 

A second interesting result following from (^) concerns the behaviour of the vortex line for small displacement d. 
By expanding E v (d/Rj_, Q, /J,) one finds that the vortex solution at d = is fully unstable if f2 < 3f2t,(/i)/5, while it 
is metastable (local minimum) 

Another important consequence of (0) is the appearance of a barrier. Even if f2 > fl v (/z) and hence if E v (d/R± , fl. u) 
is negative at d = 0, the curve ([|) exhibits a maximum at intermediate values of d between and R± (see Fig. [j]). 
The position ds and height Eb of the barrier are given by the equations 
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and 



E B = -E v (d = 0,v) 
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showing that the height of the barrier becomes smaller and smaller as f2 increases, but never disappears. Since crossing 
the barrier costs a macroscopic amount of energy, the system will never be able to overcome it and the vortex cannot 
be nucleated. 

It is finally interesting to rewrite the energy of the vortical configuration as a function of angular momentum rather 
than of the vortex displacement. This yields the expression 



E v (L z ,n,n) = Nh 



r \3/S j 



Nh) 



Nh 
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Equation ( fl3| ) emphasizes the fact that the nucleation of the vortex is associated with an increase of angular momentum 
from zero (no vortex) to Nh (one centered vortex), accompanied by an initial energy increase (barrier) and a subsequent 
monotonous energy decrease. In this form the TF-result can be usefully compared with alternative approaches based 
on microscopic calculations of the vortex energy. 



III. ROLE OF QUADRUPOLE DEFORMATIONS 



In the previous section we have shown that in order to nucleate the vortex the system has to overcome a barrier 
associated with a macroscopic energy cost. In the following sections we will show that the barrier disappears at 
sufficiently high angular velocities if we allow the system to take a quadrupolar deformation. The physical mechanism 
is especially clear for an axisymmetric trap. Above a given angular velocity, fixed by ([!]) with £ = 2, the system becomes 
energetically unstable against the excitation of quadrupole oscillations. This instability causes a continuous symmetry 
breaking which can favor the occurrence of alternative paths for the nucleation of the vortex. The appearance of this 
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intermediate quadrupole deformations is now strongly supported by experimental evidence [|3[|5J,|8J and has also been 
observed in simulations of the time-dependent Gross-Pitaevskii equation []l3|| . 

In order to describe properly the effects of the quadrupole deformation we introduce, in addition to the vortical 
field (m), an irrotational quadrupolar velocity field given by 



v Q = aV(xy) . 



(14) 



where a is a parameter. Note that vq is the velocity field in the laboratory frame expressed in terms of the coordinates 
of the rotating frame. The form (|l4|) is suggested by the quadrupolar class of irrotational solutions exhibited by the 
time-dependent Gross-Pitaevskii equation in the rotating frame J|(J • In Ref. ]2(J it has been shown that these solutions 
are associated with a quadrupole deformation of the density described by the deformation parameter 



S = 
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The parameters 5 and a characterize the quadrupole degrees of freedom that we are including in our picture. In 
the following, in order to provide a simplified description, we will fix a relationship between these two parameters by 
requiring that the quadrupole velocity field satisfies the condition: 



V • [n(r) (v e -Ox r)] = 



(16) 



where H = Clz. As a consequence of the equation of continuity, this condition implies that in the rotating frame the 
density of the gas is stationary except for the motion of the vortex core which however involves the change of the 
density at small length scales. Eq.([l6|) yields the relationship a = —Q5 which selects a natural class of paths 
that will be considered in the present investigation. In the presence of the quadrupole velocity field (0), the energy 
of the condensate in the rotating frame can then be expressed only in terms of the deformation parameter 5. Using 
the formalism of [Eol one finds the expression: 



E Q (S,Sl,e,n) = N/i 



2 l-eS-n 2 5 



2r2 



7 vi^VT 3 ^ 



(17) 



Here, e describes the deformation of the trapping potential in the x-y plane 



Vext(r) 



[(1 + e) uj]_x 2 + (1 - e) ujW + oj 2 z z 21 



(18) 



where lj± is an average transverse oscillator frequency 



(19) 



and Cl — Q/luj_ is the angular velocity of the trap expressed in units of u>±. Note that it is crucial to distinguish 
between the transverse trap deformation e and the quadrupole deformation 5 of the condensate. In Eq. (|l7| ) we have 
neglected the change of the central density caused by the velocity field ( |l4| ) . This assumption will be used throughout 
the paper pjj . 

It is useful to expand Eq. ( |l7| ) as a function of 5 in the case of axisymmetric trapping (e — 0). One finds the result 



E Q (S,n,E = 0,fj,)^Nfx 



+ s 2 



(1-2JT) +0{8 



(20) 



which explicitly shows that for ft > uj±/\i / 2 the symmetric configuration (S = 0) is energetically unstable against the 
occurrence of quadrupole deformations. As we will see in the next section this instability is crucial for the nucleation 
of the vortex. 

Our aim now is to calculate the energy in the rotating frame combining the effects of the vortex and of the quadrupole 
deformation. The parameter 6 thereby emerges as a natural degree of freedom that can be used to optimize the path 
towards vortex nucleation. In the simultaneous presence of the vortex and of quadrupole deformation the energy of 
the system is not however simply given by the sum of Eqs. Q and (p"7|). In fact, on the one hand the kinetic energy 
contains the additional crossed term mfdr [v vort cx ■ vq] n(r). On the other hand we have to take into account that 
the vortex energy (Q) and the angular momentum (^) were calculated for an axisymmetric trap and in the absence of 
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quadrupole deformation. In the presence of a deformed trap the vortex energy in the laboratory frame (Q) is simply 
generalized to 



E v (d/R x ,e,fj,) = E v (d = 0, e, ^) 
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where R x is the TF-radius of the cloud along x and we have explicitly positioned the vortex on the x-axis 122] . It is 
important to note that the first equality in (||) also holds for e ^ 0. In this case the energy E v (d — 0,e,fi) can be 
written in the form 
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E v (d = 0,e,fji) = Nhuj ± - — -y/l -e 2 log ( 1.342 
4 \x 
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which points out its dependence on the trap deformation. 

Finally, it is possible to also evaluate in a simple way the sum of the crossed term m Jdr [v vor t cx ■ vq] n(r) entering 
the kinetic energy and of the angular momentum term — m Q ■ Jdr [r x v vor t ex ] fi(r) due to the vortex. One finds 



E 3 (d/R x ,S, Q)=mjdr v vortcx • [vg - ft x r] n(r) = -Nhuj ± fl ^l-5 2 1 - (-^j 
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(23) 



which is consistent with (g) in the case of a symmetric condensate (6 = 0). In deriving result (|23j) we have used the 
condition ( |l6| ) for the quadrupole velocity field and we have integrated by parts using the expression (|3|) for v vort0 x- 
Note that the irrotational component VS of v vortex does not contribute to (E3T). 



IV. CRITICAL ANGULAR VELOCITY FOR VORTEX NUCLEATION 



The total energy of a quadrupolar deformed condensate in the presence of a vortex is given by the sum of Eqs. 
©, © and @: 

E tot (d/R x ,6,n,e,fj,) = E v (d/R x ,s, M ) + E Q (8,n,e, m) + E 3 (d/R x ,8,0), (24) 

where we have explicitly indicated the dependence of the three contributions on the various physical parameters. 
The values of the angular velocity O, the trap anisotropy e, the average transverse oscillator frequency lu±, and the 
chemical potential fi are fixed by the experimental conditions. Hence the degrees of freedom of the system with 
which one can play in order to identify the optimal path for vortex nucleation are the vortex displacement d and 
the condensate deformation 5. Note that each of the three energy contributions has a different dependence on the 
chemical potential resulting in a non-trivial dependence of the critical angular velocity on the relevant parameters of 
the system (see Eq. ( p6| ) below). Result (24) generalizes the one given in Ref. Il7|], which holds only in the limit of 



small Q where 8 ~ e. 

We assume that the system is initially in the state d/R x = 1,5 = 0. This assumption adequately describes an 
experiment in which Q and e are switched on suddenly. In this case, the condensate is initially axisymmetric (5 = 0) 
and vortex-free (d/R x = 1). Of course this configuration is not stationary and will evolve in time. In the following 
we will make use of energetic considerations in order to explore the possible paths followed by the system towards 
the nucleation of the vortex line. These paths should be associated with a monotonous decrease of the energy. In 
Figs. |2| and || we have plotted the energy surface E to t for two different values of the angular velocity Q and fixed 
e and fj,/Hu>±. In both cases the angular velocity was chosen high enough to make the vortex state a global energy 
minimum. This minimum is surrounded by an energy ridge which, for 5 = 0, forms a barrier between the initial state 
(d/Rx = 1) and the vortex state (d = 0), as discussed in section^ Moreover, the energy ridge exhibits a saddle point 
at non-zero deformation 8. The height of this saddle depends on f2, s, and fi/Hu!±. In Fig. |^ the energy at the saddle 
point is higher than the energy of the initial state and the ridge can not be surpassed. However, at higher angular 
velocities fl the situation changes. In Fig. || the saddle lies lower than the initial state. Hence in this case the system 
can bypass the barrier by crossing the saddle. The corresponding path is always associated with the occurrence of a 
strong intermediate deformation of the condensate. 

The critical angular velocity for the nucleation of vortices naturally emerges as the angular velocity O c at which 
the energy on the saddle point ((d/R x ) sp , S sp ) is the same as the energy of the initial state (d/R x =1,8 = 0): 

E((d/R x ) sp ,8 sp , Q c , s, n) = E(d/R x = 1,5 = 0, fl c , e, ll) . (25) 
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It is worth mentioning that crossing the saddle point is not the only possibility for the system to lower its energy. In 
fact Figs. H and || show the existence of stationary deformed vortex- free states which can be reached starting from the 
initial state. These are the states predicted in [M an d experimentally studied in ||,|| through an adiabatic increase 
of either Q or e instead of doing a rapid switch-on. The energy ridge separates these configurations from the vortex 
state. Still, under certain conditions this stationary vortex-free state becomes dynamically unstable and a vortex can 
be nucleated starting out from it pMITlll- The study of this type of vortex nucleation is beyond the scope of the 
present paper. 

The actual value of the critical angular velocity ft c for vortex nucleation depends on the parameters e and n/hui±. 
Fig. ^ shows the dependence of ft c on e for different choices of To lowest order in e and hujj_/(i the dependence 

is given by 



ft c (e)- 



1 1 

V2 ~ V2 
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where 



log 1.342 
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and A = 2 _5 / 4 21-\/3. This formula shows that the relevant parameters of the expansion are r/ 1 ^ 2 and e/77 1 / 4 . Fig. ^ 
demonstrates that it is applicable also at rather small values of the chemical potential. 

At e — the vortex nucleation, according to the present scenario, is possible only at angular velocities slightly 
higher than the value lu±/^/2. For non- vanishing e the preferable configuration will be always deformed even for small 
values of ft where 8 depends linearly on e. At higher ft the condensate gains energy by increasing its deformation in 
a nonlinear way (see p0[ ). Eq. (^6|) and Fig. [| shows that for non- vanishing e the saddle point on the energy ridge 
can be surpassed at angular velocities smaller than uj±/\^2. 

The above scenario seems to be in reasonable agreement with experiments. In particular, evidence for critical 
angular velocities smaller than uj±/\/2 occurring when e ^ has emerged from the experiments reported in 
A more systematic comparison with the expected dependence of ft c on the trap deformation e and on the chemical 
potential \i would be crucial in order to assess the validity of the model Q . 



V. STABILITY OF A VORTEX-CONFIGURATION AGAINST QUADRUPOLE DEFORMATION 

Once the energy barrier is bypassed, the vortex moves to the center of the trap (d/R x — 0) where the energy has 
a minimum. In an axisymmetric trap [e = 0) this configuration will be in general stable against the formation of 
quadrupole deformations of the condensate unless the angular velocity Q of the trap becomes too large. The criterion 
for instability is easily obtained by studying the ^-dependence of the energy of the system in the presence of a single 
quantized vortex located at d/R x = 0. Considering the total energy ( pij ) we find 

E to t(d/R x =0,S,n,e = 0, M ) E tot {d/R x = 0, 5 = 0, ft, e = 0, ft) + S 2 Nfi Q(l - 2ft 2 ) + + 0(5 3 ) . (27) 

Comparing Eqs. ( p7j ) with the analog expression ( p0| ) holding in the absence of the vortex line one observes that in 
the presence of the vortex the instability against quadrupole deformation occurs at a higher angular velocity given by 

n - ui b + i— )■ (28) 

If the angular velocity is smaller than (Pq) the vortex is stable in the axisymmetric configuration while at higher 
angular velocities the system prefers to deform, giving rise to new stationary configurations. The critical angular 
velocity (|28| ) can also be obtained by applying the Landau criterion ( fl|) t o the quadrupole collective frequencies in the 
presence of a quantized vortex. These frequencies were calculated in ]24| using a sum rule approach. For the m = ±2 
quadrupole frequencies the result reads 

cj ±2 =w ± V2±^, (29) 



G 



where 



is the frequency splitting between the two modes. Applying the condition ([l]) to the m = +2 mode one can immediately 
reproduce result ( p8| ) for the onset of the quadrupole instability in the presence of the quantized vortex. 

In an anisotropic trap (e ^ 0), the stable vortex state will generally be associated with a non-zero deformation 5 
of the condensate. It is interesting to note that 6 increases with the angular velocity of the trap and easily exceeds e 
(see Figs. ^ and|||). This behaviour is analogous to the properties of the deformed stationary states in the absence of 
vortices |20fl. 



VI. FINAL COMMENTS AND CONCLUSIONS 



In this paper we have developed a semi-analytic model to describe the nucleation of a quantized vortex induced 
by the sudden switch-on of the trap deformation and rotation. The trap was assumed to be of quadrupolar shape. 
The model includes the distance of the vortex line from the principal axis and the quadrupole deformation of the 
condensate as the key degrees of freedom of the system. We have shown that, for sufficiently high angular velocities, 
the instability exhibited by a condensate with respect to surface quadrupole deformations gives rise to a path for 
vortex nucleation. The energy diagram is characterized by the occurrence of a saddle point whose height can favor or 
inhibit the nucleation, depending on the value of the angular velocity, the deformation of the trap and the chemical 
potential. The scenario for the nucleation emerging from the present approach agrees with the experimental results 
recently obtained in as well as with the predictions based on the numerical solution of the time-dependent 

Gross-Pitaevskii equation ||13[. In particular, the nucleation is always associated with an intermediate configuration 
exhibiting a significant quadrupole deformation of the condensate, even in the presence of a very small deformation 
of the trap. Furthermore the nucleation is favored by increasing the values of the trap deformation and the value of 
the chemical potential. 

The energy diagram employed to calculate the critical angular velocity was based on the Thomas-Fermi approxima- 
tion to the vortex energy. This approximation is expected to become worse and worse as the vortex line approaches 
the surface region |^,^8|. Since the main mechanism of nucleation takes place near the surface, we expect that 
the quantitative predictions concerning the dependence of the critical angular velocity on the trap deformation and 
on the chemical potential would be improved by using a microscopic evaluation of the vortex energy, beyond the 
Thomas-Fermi approximation. This should include, in particular, the density dependence of the size of the vortex 
core, as well as the proper inclusion of image vortex effects. Work in this direction is in progress. 

Useful discussions with Y. Castin, F. Dalfovo, J. Dalibard, D. Feder, E. Hodby, O. Marago and A. Recati are 
acknowledged. 
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FIG. 1. Vortex excitation energy in the rotating frame (B) for an axisymmetric configuration as a function of the reduced 
vortex displacement d/R± from the center. The curves refer to two different choices for the angular velocity of the trap: 
Q, — Q v {fi) (solid line), and = 3f2 l ,(/i)/2 (dotted line), where Q v (fi) — E v (d = Q,^i)/Nh. The initial vortex- free state 
corresponds to d/R± — 1. For Q > Q v (fi), the state with a vortex at the center (d/R± = 0) is preferable. However, in 
this configuration the nucleation of the vortex is inhibited by a barrier separating the vortex- free state (d/R± = 1) from the 
energetically favored vortex state (d/R± — 0). 
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FIG. 2. Below the critical angular velocity of vortex nucleation: the plot shows the dependence of the total energy J24| ) 
minus the energy of the initial non-deformed vortex-free state E(d/R x — 1,8 = 0) on the quadrupolar shape deformation 8 
and on the vortex displacement d/R x from the center. Energy is given in units of Nhu>±. The dashed line corresponds to 
Etot — E(d/R x = 1,6 = 0) = 0, while the solid curve refers to E to t — E(d/R x = 1,8 = 0) = 0.0l5Nhuj±. This plot has been 
obtained by setting e = 0.04, n = 10 hui± and f2 = 0.64o;x. The initial state is indicated with ©, while (g> corresponds to the 
energetically preferable centered vortex state. The barrier inhibits vortex nucleation in a non-deforming condensate (5 = 0). 
The saddle point lies lower than the barrier 0. However, at the chosen Q the energy on the saddle is still higher than the 
one of the initial state 0. Note that the preferable vortex state is associated with a shape deformation 8 > e (see section |y]). 
Note also the existence of a favorable deformed and vortex-free state labeled by y Eo). 
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FIG. 3. Above the critical angular velocity of vortex nucleation: the plot shows the dependence of the total energy tej ) 
minus the energy of the initial non-deformed vortex-free state E(d/R x = 1,5 = 0) on the quadrupolar shape deformation 8 
and on the vortex displacement d/R x from the center. Energy is given in units of NfkJ±. The dashed line corresponds to 
Etot - E(d/R x = 1, 8 = 0) = 0, while the solid curve refers to E tot - E(d/R x = 1,5 = 0)= 0.015iV7w±. In this plot e = 0.04, 

= lOftt^x, as in Fig. m, and £7 = 0.7lo±. Important states are indicated as in Fig. bl At the chosen Q, the saddle lies 
lower than the initial state © allowing the system to bypass the barrier by taking a quadrupolar deformation 8 and reach 
the preferable vortex state 0. Note also the existence of a favorable deformed and vortex-free state labeled by V 
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FIG. 4. Critical angular velocity of vortex nucleation in units of u>± as a function of the trap deformation e. The long 
dashed and short dashed curves correspond to the numerical calculation satisfying condition for [i = 9 hui± and fi = 18 hu>± 
respectively, while the solid line is the analytic prediction (|2(]) evaluated with /i = 9huj±. The arrow indicates the angular 
velocity at which the quadrupole surface mode becomes unstable in the case e — 0. The value fi — 9 huj± is close to the 
experimental setting of EJ. 
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